Second-quantized Landau-Zener theory for dynamical instabilities 
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State engineering in nonlinear quantum dynamics sometimes may demand driving the system 
through a sequence of dynamically unstable intermediate states. This very general scenario is es- 
pecially relevant to dilute Bose-Einstein condensates, for which ambitious control schemes have 
been based on the powerful Gross-Pitaevskii mean field theory. Since this theory breaks down on 
logarithmically short time scales in the presence of dynamical instabilities, an interval of insta- 
bilities introduces quantum corrections, which may possibly derail a control scheme. To provide 
a widely applicable theory for such quantum corrections, this paper solves a general problem of 
time- dependent quantum mechanical dynamical instability, by modelling it as a second-quantized 
analogue of a Landau-Zener avoided crossing: a 'twisted crossing'. 



Quantum dynamical instabilities (complex excitation 
frequencies) have been associated with formation of soli- 
tons and vortices in driven Bose-Einstein condensates. 
For example, in Refs. and [||, the Bogoliubov spectra 
are examined for two-parameter families of mean field 
states. Although the physics involved is significantly dif- 
ferent (one-dimensional modulated current versus rotat- 
ing harmonic trap) , figures in both these works show dy- 
namical instabilities in narrow, finger-like regions of pa- 
rameter space (see Figure 1). We will show below that 
this is generic for dynamically unstable quasiparticle ex- 
citations in weakly interacting many-body systems. 

Such systems, especially condensates, may be suffi- 
ciently weakly damped that state engineering under ener- 
getic instability becomes feasible. Interesting structures 
such as solitons and sonic event horizons will require this 
capability. But to access the entire volume of energet- 
ically unstable control space, it will also be a common 
task to drive systems briefly through narrow regions of 
dynamical instability. Mean field theory and adiabatic 
approximations, which are important tools of state engi- 
neering, will break down during these episodes This 
paper therefore exploits an analogy between dynamical 
instability in second quantization, and avoided crossing 
in two-state quantum mechanics |4j, to determine the 
quantum effects of driving a system through dynamical 
instability. An almost identical calculation has recently 
been presented to deal with the specific dynamical in- 
stability involved in dissociation of a molecular conden- 
sate ; the present paper provides an explicit derivation 
of the result in a general context. 

Quantum dynamical instability is not always recog- 
nized as a possibility, because with a Hermitian Hamilto- 
nian, the eigenfrequencies of Schrodinger evolution must 
all be real. But the eigenfrequencies of Heisenberg evo- 
lution may be complex. It is easy to see this for the 
manifestly Hermitian Hamiltonian 




FIG. 1: Grayscale plot of sum of absolute values of imagi- 
nary parts of Bogoliubov excitation frequencies for a quasi- 
one-dimensional condensate flowing with winding number 10 
around a toroidal trap with modulated potential. Axes are 
control parameters: mean condensate density increases from 
bottom to top, severity of density modulation increases from 
left to right. Compare to figures in Refs. [nl, 



where all co-efficients are real, and a± and a\_ are ordi- 
nary annihilation and creation operators {i.e. [aj,a]^] — 
djk is the only non- vanishing commutator). The eigen- 
operators of Heisenberg evolution satisfy 

th± = [b±,H]^nJ± (2) 

{6 ± ^(52 - 72)la+ + ((5 =F ~ j^)ial 



^1^2 _^2|l/4 



(3) 
(4) 



H — {lli + S) a^_^_a+ + {S — lli) a^_a- 
-I-7 l^a^aL + a+a_ 



where we assume 7 > (without loss of generality, since 
we can always redefine a+ — > — a+). These operators 
always obey [&+,&_] = 0. As long as 7 < \S\, Sl± are 
both real, and we also have 



(1) 



7<|(5|: [b±,bi]^±sgn{S) 



[6±,feU =0 



(5) 
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And we can write 



H = sgn((5) X 



1<\5\ (6) 



for i?o = <5 — cj. This means that as long as 7 < \5\, one of 
the h± is reahy a creation operator which would have a ^ 
in more conventional notation; but apart from this trivial 
issue of labelling, we have diagonalized to noninteracting 
quasiparticles of the standard (bosonic) kind. 

If 7 > however, ^± are complex, and neither 6-1- 
nor &j_ are standard annihilation operators. They obey 



1>\5\: [hi 



A] 



±i 



(7) 



so that we must say that the canonical conjugate of 
is T*^^pj E^nd no longer coincides with the Hermitian con- 
jugate. In this case we have 



H = En 



± 



7> |<5| 



(8) 



instead of (^. One can also show that if 7 > \5\ then the 
eigenvalues of H no longer involve integer multiples of 
r2±: they remain real, but are continuous. Hence in this 
case there is no way to interpret the energy eigenstates 
as having definite numbers of any kind of quasiparticle. 

If 7 > \5\, then H is unbounded below as well as above; 
this can also occur even for 7 < \5\, if \u)\ > 5'^ — 7^- 
In these cases, (Q) can be physical only as the lineariza- 
tion of a nonlinear Hamiltonian about an excited state. 
And in all cases, since the 'anomalous coupling' term pro- 
portional to 7 in does not conserve the numbers of 
particles created and annihilated by alj_ and a±, these 
must be some kind of quasiparticle, or else there must 
be a source of atoms (such as a molecular condensate 
treated as a c- number field |^). But there are many real- 
istic situations in which unboundedness of the linearized 
Hamiltonian, and non-conservation of (quasi-)particles, 
are perfectly correct. For these (0) is universal, in the 
sense that for all eigenfrequencies, whether real or com- 
plex, pairs of modes may be expressed in this form by 
using (H) to define a±. 

Two-body interactions among particles lead to quasi- 
particle Hamiltonians with anomalous couplings as in (|^) . 
Since 7 is thus constrained to be small if interactions 
are weak, dynamical instabilities tend to arise in weakly 
interacting systems not through 7 becoming large, but 
through 5 becoming even smaller. A dynamical insta- 
bility is thus a kind of second-quantized analogue to an 
avoided crossing |^]. See Figure 2. If the many-body 
Hamiltonian depends on some number D of control pa- 
rameters, then any two modes will have frequencies ex- 
actly opposite on a {D — l)-dimensional parameter sub- 
space, and so there will be a dynamical instability within 
a thin shell around that (hyper)surface. Of course, dy- 
namical instabilities can only arise if there is energetic 





FIG. 2: At left, an ordinary avoided crossing; at right, a 
'twisted crossing' with complex intermediate eigenfrequencies. 



instability, so that some negative frequency modes exist. 
But if the system is energetically unstable, then finger- 
like (or in higher dimensioned control spaces, hypershell- 
like) regions of dynamical instability are generic. 

Since much state engineering can be accomplished 
by means of time-dependent control parameters in the 
Hamiltonian, we need to assess the effects of a time de- 
pendence that brings a system through an intermediate 
interval in which one dynamical instability exists. Fo- 
cusing our attention on the pair of modes involved in the 
instability, we can always cast them in the form of (|^), 
with time-dependent 5, 7. As long as \5\ ^ 7 we can 
apply standard methods, but a connection formula is re- 
quired to give the quantum evolution over the intermedi- 
ate interval. Within this interval, we assume the widely 
applicable case that lo and 7 are essentially constant. 



and 5 changes linearly with time. Since (a^a+ 



-) 



commutes with H , we can eliminate lo in an interaction 
picture by taking a± — *■ e^'"'a±, giving us the problem 



7(a^al + a+a_) + C^t(a^a+ + aLa_) (9) 



for a constant rate C- 

By examining the Heisenberg equations of motion 
with this time-dependent Hamiltonian we can straight- 
forwardly obtain 



lim 



a± 



it) - 



a± IN 



2CH 



2i in 



(10) 



t — *+oo 

where 



lim a±{t) = e '^'^'^Si out 



2CH 



^=F OUT (11) 



m 



■lii\2Ct\+0{t-^) 



(12) 



and a±iN,ouT are time-independent annihilation opera- 
tors, whose relationship to each other unfortunately de- 
pends on the intervening evolution at finite times. (The 
factor multiplying t in the argument of the logarithm in 
( p^ ) is not actually fixed by the equation of motion, but is 
chosen for our notational convenience below.) Of course 
the point of the limits t — s- ±00 is not literally to consider 
infinite times, but only to give the asymptotic behaviour 
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which can be smoothly matched to adiabatic approxima- 
tions vahd for C^|t| ^ 7. Hence in the Heisenberg picture 
no further definition of a±jpf is required: it wiU be given 
in terms of initial time operators by solving the evolution 
problem for t <C — C~^7- So the problem is to determine 
the OUT operators in terms of the in operators. 

This problem is directly analogous to the Landau- 
Zener problem for amplitudes in a two-state system be- 
ing driven through an avoided crossing, except that the 
amplitudes are replaced by operators, and the second- 
quantized problem allows the 'anomalous' coupling which 
would be non-Hermitian in the two-state system. Be- 
cause our problem is also linear, we can solve it in a very 
similar way. Iterating the Heisenberg equations to obtain 
a second order equation for a+ with eliminated, 



[du + 2iChdt 



2iC 



+ = (13) 



then Fourier transforming from t to frequency v, we find 
an easily-solved first order differential equation in v. 



a+{t) 



(14) 



Inverse Fourier transforming back to t then gives 



a+{t) 

A At) 



, f dv , 



(15) 



There are two independent solutions Ai^2{t)i with inde- 
pendent operator-valued co-efhcients cii,2, since the v in- 
tegral may be evaluated along inequivalent contours Ci^2 
in the complex z/-plane. Since the integrand vanishes at 
the branch point 1/ = 0, and at infinity between angles 
(0, |) and (tt, ^), integration contours can only termi- 
nate at these points (because terminating anywhere else 
amounts to replacing the integrand with something that 
does vanish elsewhere). The simplest distinct contours 
are Ci running from lower- left infinity to upper-right, 
and C2 running from zero to upper-right infinity. Since 
the integrand is analytic except at upper-left and lower- 
right infinity, and on the branch cut (which we take to run 
along the negative imaginary axis) , we are free to choose 
contours that make the evaluation easier; see Fig. 3. 

A given choice of contour yields a single solution valid 
for all but for \t\ — > 00 the integral may be evaluated us- 
ing combinations of the method of steepest descents and 
integral representations of the Gamma function. (There 
is a second saddlepoint near the origin, but as \t\ ^ 00 it 




FIG. 3: Contours used in the text, for ("^t ^> 7 (solid) 
and (^"^t <C —7 (dashed), plotted in the complex plane of 
v/{C,'^\t\). Saddlepoints are marked with crosses on the real 
axis at 2sgn(t). At left, contour Ci traverses from lower-left 
to upper-right infinity; at right, contour C2 runs from to 
upper-right infinity. 



is too close to the branch point for the method of steep- 
est descents to work, and so contour C2 must run along 
the imaginary axis until it reaches a contour of steepest 
descent through the other saddlepoint.) The asymptotic 
behaviours of the full solutions Aj{t) can thus be ob- 
tained explicitly, allowing us to determine the late time 
behaviour of the solutions whose early time behaviour 
matches (pX|). 

The contributions to the contour integrals for Aj{t) 
from the diagonal portions are dominated by the saddle- 



points vt = 2C,^t -f O(^), since for 
' dv 



^e'^ = Ids 



1 + 



2CH 



(16) 



So if the contour passes through the saddlepoint we get 
a contribution to Aj (t) 



[l + OiCt)-^]AsP = 2^/¥e*Te 




In |2Ct 



(17) 



If a diagonal segment shown in Fig. 3 does not cross 
a saddlepoint, then this segment yields only an utterly 
negligible contribution ^ e"^'' * . 

For the vertical segments along the imaginary axis in 
C2, we let 1/ = — i^/t -f and integrate 



dv 



:iSt 



1 



2{Ctf 



d^e ^ ct e 



iCi 



2C' g if-' * ^ ' 



iCt 



-r 1 
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2C 



^ (IB) 
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where in the second Hne we extend the integration hmit 
to infinity by ignoring terms of order e~^^^*^ , and drop 
the exp —i ^^-i^'i from the integrand because it yields cor- 
rections of order (^t)"^. This yields the contribution 



A^M = ^e^-(*)r ( 1 

x[i + o(C<)-^] 



17 



■ In |2Ct|] 



(19) 



From the two segments at < > where Ci hugs the branch 

cut, we gain a contribution (1 — e <^ )Aim- 

Combining these contributions and recalling our def- 
inition of 9{t) in (p^), we obtain, up to corrections of 
order (C^)^^, 



A2{t) 



t < 
t > 



(20) 



2 

iCt 



e e 



(i + iS) , 



i < 



(21) 



where the suppressed 0{Ct)^^ terms can easily be com- 
puted from (n9|, but will not be needed. 

Comparing (|2C|) with ( p^ ) and (|l^), we can read off 



' 4 e 



ai 



a2 



2 '^^^7 



-a+iN 



(22) 



e "c^ 



2C 



r 



r(i 



r(i + -^) 



using the F-function identity 

|F(l + zx)p = - 



''-IN 



e «^ 



(23) 



sinh TTX 



(24) 



Then comparing these last results with (|rT|), and noting 
that the equations for a_ and a\_ are exactly similar to 
those we have examined, we finally conclude 



Equation (|25|) is our main result. It is indeed very 
similar to the Landau-Zener formula, except for the vital 
difference that the exponents are all positive rather than 
negative. Consequently, in the state which is annihilated 
by a±77v, the expected numbers of quasiparticles after 
passing through the instability will be 



IN 



± OUT°'± OUT 



IN 



= e 



(26) 



Driving through a dynamical instability in less than the 
characteristic time scale of the instability produces only 
mild excitation, but exceeding this time generates quasi- 
particles explosively. And the crossover between these 
regimes is rather abrupt. This may set important prac- 
tical limits to the precision of state engineering, because 
controlling faster than 7 may mean not being slow enough 
to be adiabatic for other, dynamically stable modes of the 
total nonlinear system. 

Note that changing —('^ in (||), so that the sys- 

tem is stable at early times rather than late, produces 
only e"' —1 



in (25) 



And while we have analyzed 
the most general type of dynamical instability, with com- 
plex frequency and two degrees of freedom involved, the 
simpler case with purely imaginary frequency, which was 
treated in B, can be obtained by setting w — > 0, and 
a± a/\/2m (|l]). Everywhere else, one must use sim- 
ply d± ^ a instead; there are also trivial changes to 
(^ and (H). We have therefore obtained a quite general 
result that will be useful for state engineering through 
dynamical instabilities. 
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